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We consider the effects of experimental imperfections on the problem of estimation-based feedback 
control of a trapped particle under continuous position measurement. These limitations violate the 
assumption that the estimator (i.e. filter) accurately models the underlying system, thus requiring 
a separate analysis of the system and filter dynamics. We quantify the parameter regimes for stable 
cooling, and show that the control scheme is robust to detector inefficiency, time delay, technical 
noise, and miscalibrated parameters. We apply these results to the specific context of a weakly 
interacting Bose-Einstein condensate (BEC). Given that this system has previously been shown to 
be less stable than a feedback-cooled BEC with strong interatomic interactions, this result shows that 
reasonable experimental imperfections do not limit the feasibility of cooling a BEC by continuous 
measurement and feedback. 
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I. INTRODUCTION 

Ultracold atomic gases and Bose-Einstein condensates 
(BEC) are one of the premier platforms for the investi- 
gation of quantum fields £Q and precision inertial mea- 
surements [2], due to their isolation from the environ- 
ment and ability to be controlled precisely using optical 
and magnetic fields. The practical sensitivity of precision 
devices is limited by the ability to produce stable, low- 
linewidth sources [3j 0]. Measurement-based feedback 
control has shown promise for improving control of quan- 
tum systems, although early experiments [5HSj and much 
theoretical work [91413] on feedback control of quantum 
systems has been applied to relatively low-dimensional 
systems. Models of Bayesian feedback control on multi- 
mode BEC systems have shown that they can be cooled 
using feedback control for readily accessible trap parame- 
ters [MHTB] . The key assumption in most Bayesian feed- 
back control simulations is that the estimate of the state 
of the system (called the filter) is an accurate represen- 
tation of the system itself. This assumption is typically 
very robust, but obviously it can break down outside var- 
ious limits. This paper quantifies the parameter regime 
for which this assumption is valid in the context of feed- 
back control of a BEC whose position is continuously 
monitored. Our results are also applicable to the cooling 
of nanomechanical resonators |17j . the localisation of a 
particle in a double- well potential [§], and, most gener- 
ally, to the linear feedback control of a quantum harmonic 
oscillator undergoing continuous position measurement. 

The first study of feedback control on ultracold gases 
showed that it could be used to narrow the linewidth 
of a continuously pumped single- mode atom laser |18j . 
This linewidth is normally limited by phase diffusion 
caused by the strong interatomic nonlinearities. Feed- 
back significantly reduces this phase diffusion, although 



the linewidth still scales with the strength of the nonlin- 
earity. Unfortunately, semiclassical models later showed 
that continuous pumping would only produce single- 
mode operation in the presence of strong nonlinearities 
[TO1 |2"U] . This suggests that alternative methods of sta- 
bilising the spatial degrees of freedom are likely to be 
very productive in producing highly coherent atom laser 
output. 

Using feedback to control the spatial degrees of free- 
dom of a trapped atom was examined by Doherty and 
Jacobs [12j . who considered a continuous position mea- 
surement of an atom with harmonic confinement. By 
assuming an initial Gaussian state for the system, and 
applying Linear-Quadratic-Gaussian (LQG) control |13l 
Sec. 6.4], they were able to calculate the optimal cooling 
scheme even in the presence of measurement backaction. 
The evolution using an arbitrary initial state was later ex- 
amined for a linear model |16j , and nonlinear models of a 
BEC using phase-contrast imaging, which gives a contin- 
uous measurement of the density profile rather than just 
the position [TU [TS] . Cooling to near the ground state 
was still possible when using this more disruptive mea- 
surement process, and the nonlinearities in fact made the 
cooling more efficient. 

All of these simulations used a Bayesian analysis where 
the best estimate of the quantum state of the system, 
called the filter, was calculated from the dynamics of the 
system and the measurement result, and an appropriate 
control was used. This explicit separation of the system 
and filter has recently been considered in the context 
of quantum state estimation for a single qubit |21j . a 
double quantum dot charge state [35] and a BEC in a 
double well potential [53]. Except in pathological cases, 
it can be shown that the filter converges to the state 
of the system conditioned on the measurement, so typi- 
cally no distinction is made between the filter and system 
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for the purposes of discussing controllability of the sys- 
tem. However, in a situation where feedback cooling is 
competing with uncharacterised heating processes, the 
timescale of this convergence is particularly important, 
as the mismatch between the system and filter will be 
widened by the heating as it is reduced by the measure- 
ment. Although filters are typically robust to corrup- 
tion of the measurement signal, time delays and even 
miscalibrations of the system, there are obviously limits 
beyond which the controlled system will no longer be sta- 
ble. When the filter is not robust to unmodelled uncer- 
tainty, control is still possible (although not guaranteed) 
using risk-sensitive filtering |24H26j . whereby the filter 
is modified to increase robustness, but with the trade- 
off that it no longer minimises the least-squares error. 
Our results show that such an approach is not required, 
as standard least-squares filtering is robust to corrupting 
classical Gaussian noise, mismatch of system and filter 
parameters and time delay of the control signal. Further- 
more, this paper quantifies the limits of the controllabil- 
ity of a BEC, and more generally a quantum oscillator, 
with respect to these experimental inevitabilities. 

The paper is structured as follows. In Sec. [TT] we in- 
troduce our model of a system-filter separation for linear 
feedback control of a single particle in an harmonic trap 
subjected to a continuous measurement of position. Ana- 
lytic results for the system stability, average steady-state 
energy and rate of convergence to steady state are pre- 
sented in Sec. IIIIl These results are used in Sec. IIVI to 
quantitatively investigate the effects of experimental im- 
perfections on the effectiveness of the control. In partic- 
ular, we focus on the effects of corrupting classical Gaus- 
sian noise, mismatch of system and filter parameters and 
time delay of the control signal. Finally, in Sec. [V] we 
consider the specific example of a trapped, weakly inter- 
acting BEC coupled to a cavity mode, and illustrate that 
a feedback-cooled BEC is robust to experimental imper- 
fections. 



VI presents the combined conclusions of 



Sec 

these results. 



II. MODEL OF SYSTEM-FILTER SEPARATION 

Semiclassical simulations showed that atomic nonlin- 
earities improved the efficiency of the cooling [15 , which 
means that the linear case is in fact the worst-case sce- 
nario. In this limit, we can model a trapped BEC as 
a single particle, as discussed in more detail in Sec. [V] 
Therefore, we consider a quantum harmonic oscillator 
(mass to, angular frequency us) undergoing a continuous 
position measurement and linear feedback control. The 
system is described by the Ito stochastic master equation 
for the conditional density operator p t |27j : 



d Pt = -i[Hg + H<! on (t),p t ]dt + a s V [x] p t dt 
+ VasVs U [x] ptdWt, 



where Hq — (p 2 + x 2 ) /2 is the oscillator Hamiltonian, 
i? c p on (f) = u p (t)x is the linear feedback control Hamil- 
tonian with control signal u p (t), as the measurement 
strength, rjs the detector efficiency, dW t the Wiener in- 
crement which satisfies dWtdWf = S(t — t')dt, and the 
superoperators are defined as 
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T>[c] p = epe 1 — (ffcp + pc'cj 
% [c] p = cp + p& — Tr { (c + c^)/5} p. 



(2) 
(3) 



These are the decoherence and innovations superopera- 
tors respectively for any arbitrary operator c. For conve- 
nience we have expressed energy and position in units of 
hu>s and \fhj (tows), respectively. 

For closed-loop feedback control, the control signal 
must be a function of the filter 7r t , which is the best- 
estimate (in the least-squares sense) of the system p t con- 
ditioned on the information obtained from the continuous 
position measurement. Ideally, the equation of motion 
for the filter will be 0. However, since filtering requires 
some a priori information about the measurement signal, 
the system being estimated, and the choice of feedback 
control, it is possible for the dynamical equation for the 
filter to differ from that describing the system. We call 
this a system-filter separation. In this paper we consider 
three distinct experimental imperfections that would re- 
sult in a system-filter separation: 

1. The measurement signal is corrupted by classical 
noise. The measurement signal for a continuous 
position measurement has the form 



dY t = 2 v /a s r]s (x) 



■dW t . 



(4) 



where (x) P — Tr{xp t }. The Gaussian noise on 
the signal is the irreducible quantum noise from 
the weak measurement, which is required in order 
to preserve the commutation relations between the 
system operators. It gets relatively smaller as the 
strength of the measurement as is increased, but 
must always be finite. However, it is always possi- 
ble for the position measurement signal to be cor- 
rupted by classical noise due to, for example, elec- 
tronic noise. If we characterise this noise as Gaus- 
sian, then the signal fed into the filter is 



dY t =dY t + V"dWf, 



(5) 



(1) 



where dW£ is the Wiener increment describing this 
classical noise and v is the strength of the classical 
noise. 

2. Filter and system parameters differ. The param- 
eters that define the filter are the measurement 
strength op, detector efficiency r\p and oscillator 
frequency ujp- We allow these to differ from the 
system parameters as,r]s and (J g. In harmonic os- 
cillator units, it is more convenient to denote the 
deviation between the filter and system oscillator 
frequency with Awf , where ujp = ws(l + Awp). 
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FIG. 1. Illustration of the feedback control loop under a 
system-filter separation. A position measurement of the sys- 
tem, pt, gives a measurement signal dY t . This signal is first 
corrupted by some classical Gaussian noise ^udWf. The re- 
sultant signal, dYt, is then used to form an estimate of the 
system, Tit- The control signal u(t), which is a function of the 
estimate iit, is fed back into the system after being delayed 
by some time t. 



Control signal is time- delayed. In a realistic exper- 
iment it will take some finite time to measure the 
system, construct the estimate, and use this esti- 
mate to feed back to the system. This means that 
the feedback experienced by the system at time t 
will be based upon the estimate of the system at 
some prior time t — t, where r is called the time 
delay. We will assume that the control signal has 
the form 



f{t) = u{t-r) = k{pr t _ r 



(6) 



where (p)^ = Tr{j57r t } is the expectation value of 
momentum as estimated by the filter, and k > 
is the feedback strength. In contrast, if the experi- 
menter is unaware that there is a time delay, then 
the filter will model the control signal as 



S(t)=u(t) = k(p)l 



(7) 



We will assume throughout this paper that the 
feedback strength k can be accurately chosen and 
implemented by the experimentalist. 

By including these considerations, the equation of motion 
for the filter is 



dirt = —i[HQ + H£ on , TTt]dt + a F T> [x] irtdt 

+ \J a FVF H [x] Tr t (dY t - 2y/a F r)F (£)*), (8) 



where H$ = (p 2 + (1 + Aoj f ) 2 x 2 ) /2 and H* on (t) = 
u w (t)x. Note that the innovations term is proportional to 
the difference between the measurement signal dY t and 
the current best-estimate of the expectation value of po- 
sition (up to constants). A diagram summarizing the 
control scheme under a system-filter separation is shown 
in Fig. |T| 



III. ANALYTIC RESULTS FOR 
SYSTEM-FILTER SEPARATION 

Since the Hamiltonian for the system contains no terms 
higher than quadratic order in position and momentum, 
an initial Gaussian state will remain Gaussian. Further- 
more, nonclassical states evolve over a short timescale to 
Gaussian states due to environmental interactions (such 
as the measurement process) [28rf3"U] . This was verified in 
the specific context of continuous measurement of atoms 
by direct integration of the full Wigner function [16] , and 
the result is that we are free to use the Gaussian approxi- 
mation. This approximation allowed Doherty and Jacobs 
to find the optimal feedback the case when the system 
and filter are identical [T^] , but we will use it to consider 
the (robust and near optimal) linear damping when the 
system and filter differ. If p t is a Gaussian state then it 
can be precisely and uniquely represented by the Wigner 
quasi-probability distribution, 



W*(x,p;t) 



expt-fOr-x^nvO-Hs-xO] 



27iVdet(V£) 



where x T = (x,p), and 



ViJt) v x p p (t) 
V£ p {t) V p %{t) 



(9) 

(10) 
(11) 



Here 
and 



and (p) P are the means, Vj? x = (x 



viz, = (p 2 y t - m? 



VP 



the variances and 



yp = 

v xp 



(xp + px) P /2 — {x) P (p) P the joint covariance. We refer to 
the collective of these latter three quantities as the vari- 
ances. The filter 7r t can also be represented with similarly 
defined W* (x , p; t) , and V£ . 

Under the Gaussian state assumption, it can be shown 
that Eqs and pi reduce to matrix differential equa- 
tions for the means and variances [T2l [16] : 



dx? = A 



K - A^VJL^j x*dt 
T - p ^ 

dWf) 



V^V?L ff Ljx^ 
2^Y^U {dW t + v 



dx p t = (AX + K^_ T )dt + 2^V P L p dW t 

V? - A p V p + V p A T p + B p - 4 Vs V p L p L T p V p 

(l + Auj F ) 2 0^ 



where 
S = 



K = 



1 

-1 


—k 



G = 











(12a) 
(12b) 

(12c) 
(12d) 

(13a) 



(13b) 
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and = G£, A p = S, = Y,[L n Ll]'E T and D p 
S[L p Lj]S T 



Equations (12c) and (12d| are decoupled from each 



other, the equations for the means, and the feedback. 
Both are examples of Riccati matrix differential equa- 
tions, and so are guaranteed to converge to a steady state 
in the limit t -> oo [H EH: 



oo) 

v; p (t -> co) 

v: P (t -> oo) 

^(*->oo) 
oo) 

^(i oo) 



1 + Awp 

2v / 2aF??f(l + v) 

\3 



(1 + A^jr) 2 

4a_F7yi?(l + i/) 
1 

2 v / 2a s ?7s 
1 



(£f - 1) 



2y / 2as?7s 
' -(6s -1 



4a s r]s 

where for convenience we have defined 



ia 2 F r) F (l 



(i + Auj F y 



C s = A /l + 4a|7 ?s . 



(14a) 
(14b) 

(14c) 
(14d) 
(14e) 
(14f) 

(15a) 
(15b) 



In contrast, the equations of motion for the means are 
delay differential equations, and so have no analytic so- 
lution. However, to first order in the time delay, we can 
we can approximate [131 pp. 300-301] 



At 



(16) 



equation: 



dx£ « K ( (I - tK) - tA„ + AtiftV^LkLI) xjdt 



(A p - 4V^rKV^L,Lj) x£dt 



(17) 



where I is the 2x2 identity matrix. Taking the ensemble 
average of Eqs (12a) and (17), it can be shown that if 



a steady state exists, and k ^ [4Ty/a F asi] F i]sV xp (t — > 
oo)] -1 , then (see Appendix |A[) 



E[x" ] = E[x£,] = 0. 



(18) 



However, what we are ultimately concerned with is the 
average steady-state energy for the system: 



-E 



(p 2 Y +(x 2 y 

E[((i)^) 2 + ((P)^) 2 ] 
(V x P x (t^oo) + VP p (t 



The terms E [({p} 1 ^) 2 ] and E U{ 



\p ^2^ 



-►oo)). (19) 
cannot be solved 



in isolation. However, using (12a) and (17) a closed set 



of ten coupled differential equations is formed by finding 
the dynamical equation for the conditional expectation 
of every pairwise combination of (x)* , (p)*, (x)% and 
(p)j. These equations can be expressed as the matrix 
differential equation 



Vf = M f v t + b t 



(20) 



This allows Eq. ( 12b I to be approximated as a differential where 



/E[((x)H 2 J\ 

E [Wt (p)t 
e l( x )t ( x )t 

e[((p);) 2 i 

E[<P>I>>£ 

Ve[((p)0 2 ]7 



b* = 



Ml + ^J 2 
^(1 + ^Kp 

PfsV^Vp, + (3 F (1 + v)kTVZ x VZ p 

f3 F (i + v ){v; p ) 2 

HP j xp xx 

PfsVZ p V x P p + p F (l + v)kT{V£ p ) 2 
Ps(V£ x f 
PsV x p x V x r p + (3 F skTV x r x V* p 
\Ps(Vn 2 + 20 FS kTVPVZ, + (3 F (1 + v)k 2 r 2 {V^ p ) 2 J 
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(22) 



For notational compactness we have denned (3 F = 
Aa F i] F , (i s = 4a s r)s, Pfs = ^aFasVFVS and 

Qf = [/W P + (1 + A"f) 2 ] (23) 

R FS = kT(3 FS V: p - 1. (24) 



We assume this feedback strength throughout Sees IV 
and|V] Note, however, that this is only the optimal feed- 
back strength when the system and filter are identical. In 
general, the optimal feedback strength k^ pt will be some 
more complicated function, or cannot be determined an- 
alytically. 



A. Average Steady-State Energy 



In the long term limit, the matrix Mt in Eq. (20) goes 
to a constant matrix M^. Provided is invertible, 
then a unique stationary solution Voo exists, given by 



Voo = -Moo 1 boo. 



(25) 



Note that this is independent of the initial conditions of 
both the filter and the system. Thus if incorrect initial 
conditions are input into the filter, then this will not 
affect the effectiveness of the control in the long time 
limit. 

In the limit where the system and filter are identical 
(i.e. v = 0, r = 0, a F = asiVF = Vs,Auj f — 0), the 
average steady-state energy is simply that computed for 
the filter in [B]: 

EL = a F i lF {2V x l(t -> oo)V: p (t y oo) k[V^(t -> TO )] : 



2krj F 



')]■ (26) 



Minimizing Eq. ( 26 ) with respect to k gives the optimal 
feedback strength, 



"-opt 



(27) 



Thus an experimenter who believes that the filter rep- 
resents the underlying system exactly would construct 
k* pt from the filter parameters a F , rj F and Awf , assume 
v = 0, and set k = k^ pt . Explicitly, the experimenter 
would use the feedback strength 



2a F y/ffp 
1 + Auj f 




-1/2 



(28) 



B. System Stability 



Even if the stationary solution (25) exists, there is no 



guarantee that the system will converge to this steady 
state. Indeed, we expect there to exist unstable regimes, 
such as for k < 0, which result in gain (and therefore in- 
definite heating) rather than damping. Heuristically, we 
expect this to occur when the filter differs from the sys- 
tem to such an extent that the control signal is based on 
a highly faulty estimate of the system, and is therefore 
ineffective. We can quantify these regimes of instabil- 
ity by considering a perturbation v t = v t — v^. Once 
the variances have attained steady state (which they are 
guaranteed to do), the equation of motion for this per- 
turbation is 



(29) 



As outlined in standard stability analysis, the stability of 
the system of equations (i.e. whether v t vanishes, grows 
indefinitely, oscillates, etc.) is determined by the eigen- 
values of the matrix M^. If the real component of all 
the eigenvalues are strictly negative, then v goes to zero, 
and the system is stable. In this case Moo is a negative- 
definite matrix, hence invertible, and so the steady state 
(25) is guaranteed to exist. However, if the real compo- 
nent of at least one eigenvalue of M^ is positive, then 
Vt will grow, and the system will never be controlled to 
a steady state. 



C. Rate of Convergence to Steady State 

Our control goal is not only to cool the oscillator as 
close to the ground state as possible, but also as fast 
as possible. The time it takes for the system energy to 
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converge to steady state depends on (a) the initial con- 
ditions, (b) the time it takes for the variances to attain 
steady state, and (c) the time it takes for v t to attain 
steady state. The effect of the initial conditions is a tran- 
sient that dies out quickly relative to the other timescales, 
and so its effect can be neglected. The convergence rate 
of the filter variances to steady state is bounded above 
by an exponential with decay rate (see Appendix |B[) 



\/2(l + Aw F ) 
l + v 

TWCf -1-Q 



where 



n = 2(1 



(1 + Auj f ) 2 (C f - 1) 



(30) 



(31) 




.2 .4 .6 .8 1 

If 





1 




2 




3 




4 




5 




6 




7 


■ 


8 
9 


■ 


10 



^0 



An identical result holds for rC ars with the replacements 
£f £s, v — > and Aujp —> 0. This is also the rate at 
which M t Moo ! which occurs before v t reaches steady 
state. Finally, the rate at which v t converges can be 
estimated by considering the eigenvalues of M^. For if 
Xi and Wj are the eigenvalues and eigenvectors of M^,, 
respectively, then 



FIG. 2. Plot of the long time convergence rate ro as a 
function of qf and t]f, assuming optimal feedback strength 
k = h pt . 



IV. EFFECTS OF EXPERIMENTAL 
IMPERFECTIONS 



V/ 



10 



(32) 



for constants Cj. Hence, for a stable system of equa- 
tions, in the long time limit (i.e. after all transients have 
damped out) |v t | exponentially decays to zero at a rate 



max Re {Aj} . 



(33) 



In practice the variances attain steady state well before 
|v t |. Therefore (33 1 usually serves as an excellent esti- 
mate of the overall timescale. 

For the case when the filter and system are identical, 
the long time convergence rate to steady state can be 
determined analytically (see Appendix [Cl) : 



''o 



= Re|fc+ \/k 2 - 4 j 



(34) 



Setting the feedback strength to fc op t [see Eq. (27 1], a 
Taylor series expansion of tq in powers of ap^fryp gives 



= V2- 



0? F T]F 



(35) 



which shows that in the limit of sufficiently small ap 
and/or rj F , the long term convergence rate is independent 
of the measurement strength and detection efficiency. 
Outside this regime, as shown in Fig. [2] ro increases with 
both increasing a F (a stronger measurement collapses to 
steady state faster), and increasing r/p (better efficiency 
improves effectiveness of control), as expected. 



We now have the tools to consider the effect of classi- 
cal noise, inaccurate filter parameters and time delay on 
the efficacy of the feedback control. When judging the 
effectiveness of the control, we will focus our analysis on 
the following three criteria: 

1. Does the system converge to a steady state? 

2. How close to the ground state is the average steady- 
state energy? 

3. What rate does the system exponentially converge 
to the steady state? 

Doherty and Jacobs [H] and Wilson et al. [TB] have 
already addressed these questions in the context where 
the system and filter are identical. Thus in our analysis, 
where we consider a separation between the system and 
filter, results will always be given with respect to this 
identical case. 



A. Effect of Classical Gaussian Noise 

Let us examine the case where the measurement signal 
fed into the filter is corrupted by classical Gaussian noise 
of strength v [see Eq. |(5|], but there is no time delay (r = 
0) and the filter parameters agree with the system (a 
1 1 i = ag, rj = r)p = rjs and Auip = 0). Furthermore, 
the feedback strength k is set by the experimenter under 
the assumption v = [see Eq. (28)]. In this scenario, 
the matrix is the same as the case when the system 



and filter are identical, except with the replacement 77 
77(1 + v) in the steady-state variances for the filter. 
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FIG. 3. Slices of (v, a, rj) parameter space showing the aver- 
age steady-state energy (top) and long-time convergence rate 
(bottom). This illustrates the effect of corrupting classical 
Gaussian noise of strength ^/v on the average steady-state 
energy and long-time convergence rate, relative to the case 
when the system and filter are identical [see Eqs ( |26| l and 
( 34 1] . Note also that r = 0, a = hf = as , r\ = 7/f = t]s and 
Alj f = 0. 



In the absence of external heating, the corruption of 
the measurement signal by classical Gaussian noise does 
not alter the stability, and so the system is guaranteed to 
converge to a steady state (see Appendix |d|) . However, 
as shown in Fig. [3j the inclusion of this classical noise de- 
grades the effectiveness of the control, both by increasing 
the average steady-state energy (relative to E^) and de- 
creasing the rate of convergence to this steady state rel- 
ative to ro- When the detection efficiency is close to one, 
the control only loses effectiveness for weaker measure- 
ment strengths. However, as rj is reduced, both E^/E^ 
and rg/r increase with increasing v, with little regard for 
the value of a. Note that there are some regimes (e.g. 
a = 1, v = 60) where decreasing r\ actually decreases 
ro/r. This does not imply, however, that lower detection 
efficiencies result in a better control, since both r and ro 
decrease with decreasing r\ (cf. Fig.[2|. 



B. Effect of Imperfect Filter Parameters 

Now we consider the case where there is a difference 
between the filter and system parameters defining the 
measurement strength, detection efficiency and oscillator 
frequency (aj? ^ as,)]F 7^ ^s, Auip =^ 0), but there is no 
time delay (r = 0) or classical Gaussian noise (y = 0). In 
this scenario, there exist some regimes where the filter's 
estimate of the system is sufficiently inaccurate that not 
only is the feedback effectiveness reduced, but is now 
entirely ineffective. In this regime, the oscillator heats 
indefinitely, and fails to converge to a steady state. Plots 



showing regions of stability, average steady-state energies 
and convergence rates relative to and ro, respectively, 
for slices of parameter space (ap, as, Tjp, rjg, Awp) are 
shown in Fig. [4] 

Before we proceed, let us consider qualitatively what 
effects a mismatch between the filter and system param- 
eters should have on the effectiveness of the control. For 
feedback to be effective, the rate at which energy is re- 
moved from the oscillator must (eventually) balance the 
rate at which the measurement backaction causes heat- 
ing. The rate at which the system heats is fixed by as- 
However, the cooling rate is strongly dependent on the 
filter. All three filter parameters are used to determine 
the feedback strength (which is optimal if the system and 
filter agree) . Imperfect selection of parameters will result 
in sub-optimal feedback, which in this context will result 
in under or over damping of the oscillator. Furthermore, 
the quality of the estimate (p) t depends on ap and r\p. 
Larger values of these parameters imply that the filter 
overestimates the momentum diffusion and information 
rate of the position measurement. 

Let us consider the case Aup — (see (a) and (b) in 
Fig. [1]). Firstly, when ap > as and/or t]f > ijs the 
system converges to a final steady-state energy approxi- 
mately equal to the identical system-filter case. Further- 
more, the convergence rate r ss ro. It seems, therefore, 
that over-estimating the measurement strength and de- 
tection efficiency has almost no effect on the effectiveness 
of the control. Secondly, if there is a weak measurement 
strength and/or low detection efficiency, then the con- 
trol is robust to imperfect guesses of as and 775. We 
have seen that filter-based estimation is robust to uncor- 
related errors with zero mean, and so it is unsurprising 
that it is also robust to an incorrectly estimated detection 
efficiency. The robustness with respect to an incorrectly 
characterised measurement strength, which causes mo- 
mentum diffusion, was less clear. 



The regime where Atup 7^ behaves somewhat dif- 
ferently [see Fig. |4jc)]. When the harmonic oscillator 
frequency is mischaracterised, then even if the filter has 
a perfect guess of the initial state, then the mean values 
will diverge from the system. The phase of the feed- 
back will then modulate. Measurement will update the 
filter to help preserve the phase, but there are obvious 
regimes of instability if the trap is mischaracterised by 
a significant fraction. There is a second effect, whereby 
the strength of the damping is chosen by the filter, so 
the feedback response is under- and over-damped on the 
lower and upper sides of the resonance Aujp = respec- 
tively. The combination of these two effects lead to an 
asymmetric response with respect to the sign of Au>f be- 
tween the filter and system. Note that there are unstable 
regions lying on both sides, and a valley of stability near 
the resonance. 
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FIG. 4. Slices in (olf, as, Vf, Vs, Aljf) parameter space 
showing the average steady-state energy (top) and long-time 
convergence rate (bottom). This illustrates the effect choos- 
ing filter parameters «f, T}f, Auf that differ from the system 
parameters has on the average steady-state energy and long- 
time convergence rate, relative to the case when the system 
and filter are identical [see Eqs ( |26[ ) and ( |34[ )]. The white 
dashed lines indicate the points where the system and filter 
are identical. Note also that r = v = 0. 
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FIG. 5. Slices in (r, a, rf) parameter space showing the aver- 
age steady-state energy (top) and long-time convergence rate 
(bottom) under short time delay approximation (16 1. This 



illustrates the effect of a time delay on the average steady- 
state energy and long-time convergence rate, relative to the 
case when the system and filter are identical [see Eqs ( 26 I and 
pi|]. Note also that v 

Aujp = 0. 



0. 



OIF 



as, V = r/p — rjs and 



C. Effect of Time Delay 

We now consider the case of nonzero time delay r 
when the filter and system parameters agree, and there 
is no classical noise (a = ap = as, rj = rjp = rjs and 
v = Awp = 0). By computing the eigenvalues of Meo, 
and using Eq. (I25|, we examined the effectiveness of the 
control (see Fig. [5J). The results show that the control is 
robust to short time delays, but very quickly loses effec- 
tiveness and becomes unstable as r increases. However, 
the boundary of instability is in fact due to a breakdown 
in the short time-delay approximation (16). We used a 



numerical solution of the stochastic equations to deter- 
mine stability for longer delay times. 



We computed the average steady -stat e ener gy ( 19 1 by 
numerically integrating equations (12a) and (12b) with 
a fixed step 4th order Runge-Kutta algorithm. This was 
done using the software package xmds2 [35]. The results 
of this numerical analysis are shown in Fig. [6] It was 
assumed that the variances had reached the steady-state 
values (14), and that the means were initially zero. For 
the parameters considered, steady state occurred some- 
where between t = 10 and t = 1000 (where t is in units of 
Wg 1 ). Instability in the system is easily recognised by ex- 
amining the energy, which when unstable grows exponen- 
tially without bound. We categorized a point as unstable 
if at the end of the interval of integration E^/E^ > 100. 
This heuristic worked well in practice, as most of the un- 
stable points sampled grew to energies 10 6 — 10 s times 
larger than E^. 

As can be seen from Fig. [6j the control is stable for 
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FIG. 6. Slices in (r, a, rf) parameter space showing the aver- 
age ste ady- state energy, computed via numerical simulation 



of Eq. ( 12a 



Points required averaging over 10 5 — 10 6 trajec- 
tories, and have a standard error no greater than 10%. Note 
also that v — 0, a = of = as, rj = r/F — gs and Aujf = 0, 
and that the r axis has a different scale to the plots shown in 
Fig. |5J 



larger regimes of parameter space than was predicted 
under the small time delay approximation (c/. Fig. [5]). 
Instability begins to occur around r > 0.6, which corre- 
sponds to the feedback lagging the oscillator by > 35°. 
For weaker measurement strengths, the feedback can lag 
the oscillator a larger amount - for some parameters 
> 60°. Furthermore, the short time-delay approxima- 



tion ( 16 ) predicted that as r increases there will always 



be smaller values of a for which the system is stable. The 
numerical analysis shown in Fig. [6] disagrees, and seems 
to indicate that after around r ~ 1.4 feedback will never 
cool the oscillator to a steady state. 



V. EXAMPLE SCENARIO: A 
FEEDBACK-COOLED BOSE-EINSTEIN 
CONDENSATE 

We now consider a scenario that incorporates all the 



experimental imperfections examined in Sec. IV simulta- 
neously; a feedback-cooled non-interacting BEC. A po- 
sition measurement could be engineered by placing the 
condensate in a cavity, probing the cavity with a laser 
off- resonant with the transition of the BEC, and measur- 
ing the output from the cavity (see Fig. [7]) . More pre- 
cisely, this gives a measurement of the centre of mass 
position X = ~Y^;Xi/N a , where Xi is the position of 
the «th atom and N a is the total number of atoms in 
the condensate. The derivation showing that this phys- 
ical situation results in a position measurement follows 
that in Doherty and Jacobs with the replacements 
x — > X and go — > %/ N a go (where go is the cavity QED 
coupling constant). In this case the system state is de- 
scribed by conditional master equation Q with x — > X, 
p — > P = '^2 i Pi/N a (for pi the ith atom's momentum) 
and measurement strength 



FIG. 7. Diagram depicting the feedback control of a BEC 
via a cavity-mediated measurement of the centre-of-mass po- 
sition. 



where fco = 2tt/X is the wave vector of the probe laser, 
n the steady state average photon number in the cav- 
ity in the absence of the atomic sample, n the cavity 
linewidth, A the detuning of the probe laser from the 
cavity resonance and Xho = ^/^/(tows). Note that this 
choice of as assumes that X has been written in units 
of iho /\/N a , which is the natural lengthscale to use in 
this situation. 

We consider a condensate of Rubidium 85 atoms, as 
this allows us to turn off the interatomic interactions 
via a Feshbach resonance (33) . A similar situation has 
been engineered in an optomechanical system, allowing 
a measurement of the centre of mass position for a (non- 
condensed) sample of cold atoms [51] . We set the system 
parameters based loosely upon those used in that exper- 
iment: N a = 10 4 , A = 780 nm, w s = 2tt x 110 kHz, 
. 9o = 2tt x 12 MHz, re = 2tt x 2 MHz, A = 2tt x 20 GHz, 
fi = 0.8 ( a s « 0.1) and t) S = 0.16. 

Finally, as a highly conservative upper bound, we will 
assume that the estimate of the system parameters is in- 
correct by 100%: a F = 0.05, n F = 0.08 and Auj f = 1. 
Furthermore, the measurement fed into the filter has clas- 
sical Gaussian noise of strength v — 10, and the control 
signal is delayed by time r = 0.1. 

Combining all of these experimental imperfections, nu- 
merical simulation of Eqs ( 12 ) shows that the control still 



as 



^N a g^nx^ 
ujskA 2 



(36) 



brings the BEC to a steady state, as is shown by the red 
curve in Fig. [8] This of course is in excellent agreement 
with the analytic results. Compared to the case where 
the system and filter are identical, the system has an av- 
erage steady-state energy roughly 4.2 times larger and 
convergence time ~ 100 times longer. 



VI. CONCLUSIONS 

In this paper we have investigated the linear feedback 
control of a quantum harmonic oscillator undergoing a 
continuous position measurement under the more realis- 
tic scenario of a filter-system separation (see Fig. [IJ . In 
particular, we considered the effectiveness of the control 
(i.e. whether the oscillator cools to steady state, and if 
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FIG. 8. Numerical simulations of the average system energy 
as a function of time for (blue) identical system and filter 
(ap — a$ = 0.1, rjs = r]F = 0.16, and Awjr = v = r = 0) 
and (red) a system-filter separation (aj = 0.05, as = 0.1, 
r) F = 0.08, rjs = 0.16, Alj f = l,v= 10, and r = 0.1). Initial 
conditions for system and filter were identical, and {x)g = 2, 
(p) P = 1, V&(0) = 2,V X %(0) = 0.25, and V&(0) = 1. Sim- 
ulations were over 10 paths, and so the standard error in 
the means is less than the line thickness on the plots. The 
black, dashed horizontal line is the average steady-state en- 
ergy for the system-filter separation (parameters as above) 
analytically computed using Eq. (25 1. 



so to what energy and at what rate) when (i) the mea- 
surement signal fed into the filter is corrupted by clas- 
sical Gaussian noise, (ii) the filter parameters govern- 
ing the measurement strength, detector efficiency, and 
oscillator trapping frequency differ to those of the sys- 
tem, and (iii) the feedback control signal is time delayed. 
Although our investigation has found regions where the 
control is clearly ineffective, these only occur for seri- 
ous mismatches between the filter and system. This was 
illustrated by considering the specific example of cool- 
ing a BEC. Overall, we conclude that, under the likely 
experimental imperfections that result in a system-filter 
separation, the control scheme is robust. 



Assuming the existence of a steady state on average, 
the ensemble average of Eqs (12a) and (17 1 gives 



<l p 1 

v ooJ ' 



= M 1 E[x^]+M 2 E[x^ o ] 
= M 3 E [x£,] + M 4 E [x£ 

where 

Mi 
M 2 

M 3 

M 4 

Now 

det (Mi) = (1 + Alo f ) 2 + Aa F r] F [kV^ x (t -> oo) 

+ (l + Au JF ) 2 v: p (t^^)}, (A3) 
det (M 4 ) = 1 - 4:kT^/a F a s r] F T]sV^ p (t -> oo). (A4) 

which are both guaranteed to be nonzero for k > and 
k 7^ [4:Ty/a F asi r i F i]sV xp (t — > oo)] -1 . We can therefore 
write 



A T + K - 4j 7F V^L 7r Lj 

K ((I - rK) - rA, + A Vf tV^ULI 
A p - V^srKVjI^L 



T 



(Ala) 
(Alb) 



(A2a) 
(A2b) 
(A2c) 
(A2d) 



and 



E[x^]--Mr 1 M 2 E[xg ] (A5) 



0= (l 2 -Mr 1 M 2 M7 1 M 3 )E[x^]. (A6) 



It is easily checked that the matrices Mi X M 2 and 
(I 2 - M7 1 M 2 M7 1 M 3 ) are nontrivial, which implies 
that E [x^J = E [xgj = 0, as required. 



Appendix B: Derivation of Bound on Convergence 
Rate for Variances 



Consider the difference between the matrix of variances 



and its steady state, V t = - V£,. Using Eq. \\2c 
and 
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= A ff V" 



■ V T A 1 



OO ' 



-^(l + i/jv^L^y 

it can be shown that 



(Bl) 



7 T 



VJ = A w \ t + V t A 7T + Arj F (l + v)V t L n L^V t , (B2) 

where A T = A, - 4r) F (l + ^V^L^L^. This is a Lya- 
punov differential equation, and has solution 



4r/ F (l + v) 



dS e A "(*- s )vX^v; e A ^*- s ). 



(B3) 
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The two eigenvalues of are 
= (1 + Aw F ) 

± V2(l + v) 
where 

Q = 2(1 + !/) [1 

The real component of these eigenvalues is always neg- 
ative, and so the integral on the right hand side of 
Eq. ( B3 1 is bounded from above by some constant matrix 
C. Hence 



± v/£ F - 1 - n (B4) 



(1 + A^ F ) 2 (£ F -1)] . (B5) 



Appendix D: Proof of System Stability When 
Measurement Signal is Corrupted by Classical 
Gaussian Noise 



In Sec. IIV Al we claimed that when there is no time 
delay (r = 0) and the filter parameters agree with the 
system (a = ap = a$, rj = rjp = rjs and Aw F = 0), 
but v 7^ 0, the system and filter always converge to a 
steady state. This can be straightfowardly shown by ex- 
amining the matrix M [see Eq. ( 22 1] . In this regime, the 
eigenvalues of are 



vr<e 



exp 



max Re { XA t 

j={+,-} 



= exp(Re{2A+}t)C, 



(B6) 



for some constant matrix C. This gives the result (30 1. 
A similar argument allows one to bou nd V t . As an aside, 
one can see immediately that bound ( B6 ) is highly insen- 



sitive to the initial condition, which can only change the 
bound by a multiplicative factor. 



Appendix C: Derivation of Convergence Rate ro 

Since we are assuming that the filter and system are 
identical, we drop the superscripts and subscripts tt, p. 
The stochastic differential equation for the means [see 
Eq. ( 12b )] is an Ornstein-Uhlenbeck process. The corre- 



Ai = 

A 2 = 

A 3 = 
A 4 = 

A 5 = 

A 6 = 

A 7 = 

A 8 = 

A 9 = 

Aio = 

where 



-k 
-k 
-k 



\fk 2 



2 -4 
Vk 2 -4 

-4ar,V* x - 2^a V (a V (V- x ) 2 - V*) - 1 (Die) 



(Dla) 

(Dlb) 

(Die) 
(Did) 



-(k + 4ar,V: x )- VC- 
-(k + 4ariV£) + y/c: 

-(k + Aar 1 V x l) + ^ 



-Aa V V: x + 2^4ar, (a V (V* x ) 2 - V&) ~ 1 (Dlf) 
1 (Dig) 

(Dlh) 
(Dli) 
(Dlj) 



lation function is [35l p. 109] 



E 



ml 



E [xf x t ] 



C± = (k + iarjVZJ 2 - 8 [1 + ar, (kV x l + 2V? p )] 



x n e 



± yj(4 - k 2 ) [1 + Aar, (V* - tw?(V*) a )] , (D2) 
and we have assumed the variances have attained their 



f. ciiiu wt iidvt: ct&suiiicu Lilt: vcllicuil.vjs iidvu dLteuiiuu Liiuii 

Aij / ds h T V T e~ AT( ' t ~ s ^e~ A< ' t ~ s ^V a L steady-state values given by Eqs ( 14 ). We will now show 
Jo that the real component of these eigenvalues is alwavs 



< x^e- A *e- At x + C, (CI) 

where C is some positive constant that will have no effect 
on the convergence time. The first term on the right hand 
side can be explicitly computed, which gives 



E 



-ki 



< 



k 2 -4 
+ k^k 2 -±({x)l-{p) 

(ty/k 



((P> 

+ k (x) (p) + (x) 



x sinh 



kU(x) Q {p) + k((x)l + (p) 2 ] 



that the real component of these eigenvalues is always 
negative, implying that the system always converges to 
the steady state. 

• Ai, A2, A4, A5: For a > 0, r\ > and v > we have 
V xx > 0,V xp > and k > 0. Under these condi- 
tions, a quick inspection shows the real component 
of these eigenvalues to be strictly negative. 

• A 3 : For k 2 < 4, A 3 = -k+u/i - k 2 => Re(A 3 ) = 
—k < 0. For k 2 > 4, \Jk 2 - 4 < k, and therefore 
Re(A 3 ) < 0. 



x cosh 



(t^k 2 



C. 



(C2) 



• Xe'- Using Eqs (|14|) it can be shown that 

av(V xx ) 2 



V"* = - 

xp 



(1 + 2j/)(1 + Aw f ) 2 (£ f -1) 



8077(1 + v) 



There are three exponentially decaying rates here: k and 
Re [k ± \Jk 2 — 4}. The slowest decaying rate is 



< 0. 



(D3) 



r = Rc|fc+v / fc 2 -4|. 



(C3) 



Therefore Aarj (ar](V* x ) 2 - V£ p ) - 1 < 0, which 
implies that Re(A 6 ) = —4m]V xx < 0. 
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A7, As, Ag, Aio: It is enough to show that 
Re (vC±) < (k + AariV^x)- First, note that we can 
rewrite the expression under the square root in £± 



a.s 



K = (4 - fc 2 ) [1 + 4m? (V: p - a V (V: x f)] 
= 4 [1 + ar, (kV x l + 2V£)] 2 -[(k + 4my^) 2 
+ia V V: p (fc 2 + Aa V (V: p + kV£))] (D4) 
<4[l + a V (kV: x + 2V: p )] 2 . (D5) 



If fc 2 < 4 then n > and so C± is real. In this case 



C- < ^(fc + 4m^) 2 - 8 [1 + ar, {kV? x + 2V%)] 
Kk + AarjV^, (D6) 



and 



^ < ^{k + 4a V V* x ) 2 - 6 [1 + ar, (kV« x + 2V£)] 
<fc + 4£» 7 V£, ) (D7) 



where we have used inequality (D5 1 to bound 



If fc 2 > 4 then k < 0, and so £± is complex. In order 
to compute \/C±, we note that if (± = a± + ib±, 
then 



Re 
Im (7C±) 



(V C: ) = J 2 \/\ +'' 2 ■ 



sgn(b±) 
72 



a+ + b± - a± . 



(D8) 
(D9) 



For the case under consideration 



fore, with some algebraic manipulation, 



= {(fc + 4myV^) 2 - 8 [1 + my + 2V^)] }' 

+ (fc + 4myVJJ 2 [1 + 2a V (kV£ x + 2V£)] - M 

< (fc + 4a V V x l) 4 + 64 [1 + ar, (kV^ + 2V? p )] 2 
- 14 (fc + 4ar ? ^ x ) 2 [l + my + 2t£)] 

< (fc + 4a V V: x ) 4 + 64 [1 + my (fcV£ + 2V^)] 2 
+ 16 (fc + 4myV^) 2 [l + ar, (fcV£ + 2^)] 

= {(fc + Aa V V: x f + 8 [1 + ar, (fcV^ + 2V£)] }' 

(D14) 



Further algebraic manipulation of inequality ( |D14 ) 
yields 



\ja 2 ± +b 2 ± + a± < 2 (fc + 4myV7J 2 



(D15) 



and therefore Re(y^±) < (fc + Aar,V xx ), as re- 
quired. 



a ± = (k + Aar,V: x f - 8 [l + my + 2^,)] (D10) 

b± = ±y/=K. (Dll) 

We can write 

b 2 ± = (fc + 4my^) 2 [l + 2ar, (kV x l + 2V£ p )} - /x, 

(D12) 

where 

/i = 2 (fc 2 - 2) + 2my [4 (fc 2 - 2) F/ p + fc (fc 2 + 4) V£\ 
+ 4a 2 r, 2 V« x [8kV: p + (8 + 3fc 2 ) l£] 
+ 32aV(^) 2 + 2 ^p) • (D13) 



is guaranteed to be positive since fc 2 > 4. There- 
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